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We study a class of FRW spacetimes with a nonminimally coupled light massive scalar field. 
Values of the coupling parameter ^ < 0 enhance long range power in the vacuum expectation 
value of the energy-momentum tensor (T);,y) and fundamentally alter the nature of infla¬ 
tionary particle production: the energy density of created particles behaves like an effective 
cosmological constant, leading generically to Ulm < 1 in clustered matter and providing a 
possible resolution of the “fl problem” for low density cosmological models. 
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In this Letter we discuss particle creation of light non¬ 
minimally coupled scalar fields due to the changing geom¬ 
etry of a spacetime which underwent an early inflationary 
phase. Nonminimally coupled fields have arisen in di¬ 
verse cosmological contexts, e.^, density perturbations 
from strongly coupled scalars |]y, the possibility that a 
very light scalar particle such as the axion may couple 
nonminimally to gravity ^ and production of primor¬ 
dial magnetic fields from nonminimal electromagnetism 
[^. Ultra-light scalars (with m ^ H) have been discussed 
in the context of pseudo-Nambu-Goldstone bosons |^. 

As shown below, for negative values of the coupling 
the dominant contribution to both ($^) and {T^if) comes 
from modes having wavelengths larger than the Hub¬ 
ble radius which makes both quantities formally infrared 
(IR) divergent. IR divergences cannot be renormalized 
away, instead a physically motivated IR cutoff has to be 
invoked |^. The IR-regularized {T^^) describing particle 
creation is finite, and behaves like an effective cosmo¬ 
logical constant as the Universe expands, consequently 
the energy density of created particles can dominate the 
matter density leading to Ulm < 1 in a flat Universe. 

We consider a spatially flat FRW model with expansion 
factor either de Sitter-like a oc e^*, or power law a ocP’. 
Massive free scalar fields satisfy the wave equation 

[□+ = 0 ( 1 ) 

where R is the Ricci scalar and ^ parametrizes the cou¬ 
pling to gravity, ^ = 0, 1/6 corresponding to min¬ 
imal and conformal coupling respectively. In a spa¬ 
tially flat FRW Universe the field variables separate and 
= (27r)“^/^(/)fc(77) The comoving wavenumber 

k = 2'Ka/X where A is the physical wavelength of scalar 
field quanta. Defining the conformal field Xk = and 
using R = 6d/a^ (differentiation beiM with respect to 
the conformal time 0 = 1 dt/a), Eq. (Q) leads to, 

Xfe + - (1 - 6^)d/a]xfe = 0 . (2) 

For de Sitter space (|^) has the exact solution 

Xkiv) = + C 2 ^H‘d\kri) , ( 3 ) 


where p?' = 9/4— 12^ —ci = C 2 = 0 gives 

the state associated with the Bunch-Davies vacuum . 

For ultra-light fields m/iJ ~ 1 today, corresponds to 
m/iJ ~ 10“®° during inflation, hence the field modes can 
be treated as being effectively massless. Setting m = 0 
in (i and specializing to a power law expansion 


a = (t/to)^ = 



(l-2i/)/2 


(4) 


where 2v = (1 —3p)/(l —p), (the inflationary rangep > 1 
corresponds to v > 3/2 and for a oc e^‘, v = 3/2), 
we find exact solutions of (§) having the form (^ with 

_ 1/4 = _ i/4)(-i _ Pqj. ^ ^ i/g^ ^ ^ 1/2 

while for ^ = 0, /r = v. The choices ci = C 2 = 0 

correspond to the adiabatic vacuum state. 

To study quantum fluctuations we define the operator 


4>(a:) 



afc$fc(x, rj) +al<I>l{x,r]) 


(5) 


where Ofc, are annihilation and creation operators 
[ak,a\,] = Skk', defining the vacuum state afc|0) = 0 Vfc. 
The two-point function 

($(a:)$(a:'))«ac = J d^kC^'^^-^'^cfkioWkiv') (6) 

is IR divergent over a certain range of p values. To see 
this, substitute Eq. (^) in Eq. (^) , using (fk = Xk/ a and 
the small-argument limit of the Hankel functions 




{kri/2Y 


fcr;<l r(l -b p) TT 


± Air) I Y 




(7) 


The integral controlling the presence of IR divergences is 


J dkk^k C2p . (8) 

For the adiabatic vacuum state d), IR divergences arise 
when fjY > 9/4. For ^ = 0, IR divergences are encoun¬ 
tered for p > 2/3 {p Y 1) The situation is substan¬ 
tially different for ^ 0 (Fig. 1), IR divergences being 
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present over a wide range of expansion rates: p < 0 ; 
p > 1/2 (p 7 ^ 1); 0 < p < 1/2. Since there is no particle 
production when p = 0 , 1 / 2 , ^ = 1 / 6 , these special cases 
are free of IR divergences. 



FIG. 1. IR divergent regions (shaded) of (3>^): The special 
casesp = 0, 1/2, 1 and 1/6 3/2 have no IR divergences. 


Significantly, {T^^) can also be IR divergent: The 
(bare) energy density in a spatially flat Universe is 


P= (Too) = 






dkk^ 




dkk^ 


H i/>: 


2 I , M*k+^k<l)*k 


+ 


47r^ 


dkk^\c^k? . 


(9) 


Restricting attention momentarily to ^ = 0, m = 0, we 
conclude that (Tqo) is IR divergent for \v\ > 5/2, *.e., 
;U4<p<2 (pT^I). [We correct a minor error in Ref. 

who quote 2/3 < p < 2 (p 7 ^ 1).] This is a much 
smaller range than for the two-point function which is 
divergent for p > 2/3 (p 7 I 1). In contrast, the key as¬ 
pect of nonminimal coupling is that (Too) contains terms 
proportional to ^(4>^), and consequently (for ultra-light 
fields) (Too) is IR divergent over the same range of pa¬ 
rameters as ($^). 

The curing of IR divergences requires that mode func¬ 
tions be modified in the IR limit. The only freedom 
to accomplish this in Eq. (||) is to change the behav¬ 
ior of |ci — C 2 I as k —> 0: There are no IR diver¬ 
gences if limfe^o |ci — C 2 p oc ^ 0 (main¬ 

taining C 2 = 0, Cl = I at large k). One way to de¬ 
termine Cl and C 2 is to assume the existence of a pre¬ 
inflationary radiation-dominated phase ^, 0 , as a result 
|§, |ci — C 2 p — [1 -f 47 r^C(fc 77 o )^~^'^']~^5 where 770 marks 
the onset of inflation. Finiteness in the ultraviolet is 
achieved by imposing a cutoff at the horizon scale. It 
then follows that 


($ 2 ) = 





( 10 ) 


where C = Cpg The behavior of (d*^) depends cru¬ 
cially upon the value oi v — \fi\. For <^ = 01 — 


C dk which gives for exponential inflation 

(iz = 3/2) the standard result (d)^) ~ At/( 47 r 2 ) |^. In 

the case of power law inflation v > 3/2, and (d)^) freezes 
to a large value at late times ||j^. With negative val¬ 
ues of p > V, (4)2) grows with time approaching the 
asymptotic form: 


($ 2 ) 


C 

17^0 2|p| — 3 


( ^ 
V V 






( 11 ) 


at late times. For ^ —1/6 and » 1/4 we have 

IpI « 7 )y^ 6 |^|, which substituted in ( 0 ) gives 

(4>^) oca^ c= 1 ) > 1 ( 12 ) 


Thus ^ < 0 can greatly accelerate the growth of fluctua¬ 
tions in inflationary models. Even for minimal coupling, 
^ = 0, IR finite physical quantities, e.g., {^{x)^{x')) 
and (Too), remain sensitive to the presence of long range 
power in the field modes. This is reflected in the growth 
of scalar fluctuations, generation of density fluctuations, 
and the quantum creation of gravitational waves. 

We now examine particle production in a Universe that 
inflates and then transits to a matter dominated regime 
of expansion. To do this, we return to Eq. (j^): this equa¬ 
tion closely resembles the one dimensional Schrodinger 
equation, the role of the “potential barrier” V(x) being 
played by V(p) = —m^a^-l-(l — 6^)d/a jl^. The form of 
the barrier is shown in Eig. 2 assuming ^<1/6, to~ 0. 
The process of super-adiabatic amplification of zero-point 
fluctuations (particle production) can be qualitatively de¬ 
scribed as follows: the amplitude of modes having wave¬ 
lengths smaller than the Hubble radius decreases confor¬ 
mally with the expansion of the Universe, whereas that 
of larger-than Hubble radius modes freezes (if ^ = 0) 
or grows with time (^ < 0). Consequently, modes with 
^ < 0 have their amplitude super-adiabatically amplified 
on re-entering the Hubble radius after inflation (Fig. 2). 

The non-vacuum state of the scalar field at late times 
(77 > I 770 I) is described by a linear superposition of posi¬ 
tive and negative frequency solutions 

f>out(k,v) = -h (13) 

where 0+" = (,/W^/2)(p/po)‘^Ill^^j^\kp). [a oc tP, p < 
1; 2P = (1 - 3p)/(l - p), -3/2 < P < -1/2 ^ 
2/3 > p > 1/2 corresponding to matter with an equa¬ 
tion of state 0 < Tu < 1/3, ( 2 ^ —1/4 = (P^ — 1/4)(1 — 6 ^).] 
The transition from inflation to a radiation/matter dom¬ 
inated epoch is marked by 770 , iT/^ = I /770 being the 
Hubble parameter at reheating. [We assume < |^|i? 
or m?/H'^ < 61 ^ 1(722 —1/4)/(j/—1/2)2 allowing us to treat 
field modes as being effectively massless.) 

The Bogoliubov coefficients a and /3 are deter¬ 
mined by matching (/)//(fc, 77 ), (f>f^{k,p) given in and 
4>out{k,p), ^out{k,p) at 77 = 770 . For modes with kpo < 1 
we obtain 
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a + p = A-T{fi)T{l + fl) 

TT 


kyo 
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-(m+m) 


+B 

a-P = C 
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r(i + M) 

in 


r(M)r(i + Ai) 
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/i+/i 


+13 


r(+) kyo 

m 


lap- 


= 1 


(14) 


where, A = {p,+ \D\ + v-y,)/2p,, B = {ji+\D\ + u + y)/2ji, 
C = (y + y, + \v\ — y)/2y, D = {y — \v\ + y — v)j2y. 
(Indices y refer to the inflationary epoch; D, y to the 
matter/radiation dominated epoch.) 



FIG. 2. The time-like “potential barrier” V{rj) for infla¬ 
tion followed by a matter dominated epoch. 


Since |/3p oc and y = y{v,^), the num¬ 

ber density of created particles is sensitive to: 1) the 
inflationary expansion rate parametrized by 2) the 
equation of state after inflation, parametrized by v, 3) 
the coupling to gravity More particles are created for 
^ < 0 than for ^ = 0 (see Fig. 2 and Ref. |^). 

From (|^) and (14) we find that on larger than horizon 
scales {ky < 1), 

(j)out ^ i A^/yo/4:TTr{y){kyo/2)~>^{y/yo)>^~^''^ (15) 


from which follows the important observation that 
{^^)out = l/27r^/ dk\(j)out\‘^k^ has exactly the same in¬ 
frared properties as (<l>^)i„ and that on scales larger 
the horizon ($^)out grows with time at a rate determined 
hy y — \v\ >0 {\D\ = 3/2 in a matter dominated Uni¬ 
verse). This result bears directly on the energy density 
of created particles (Too)out, which can be determined 
from (^) after substituting cj) (pout ■ Eq- (§) informs us 
that the main contribution to {ToQ)out comes from long 
wavelength modes and, for ^ < 0, {TQo)out is IR divergent 


if the field is effectively massless. A radiation dominated 
phase prior to inflation leads to an effective IR cutoff 
kmin — yo~^ in (Too), since IR divergent states cannot 
arise from IR finite initial conditions [|l4| . In addition, a 
high frequency cutoff appears due to suppression of par¬ 
ticle creation at large k [Q. For particles created during 
inflation, this cutoff is set by the Hubble parameter at 
the end of reheating and the commencement of radiation 
domination: kmax — Hrh — yo^ 

The integration limits in are therefore For 

m/H < 1, (Too) is dominated by modes larger than the 
Hubble radius, thus the integration limits are effectively 

leading, for a matter dominated Universe, to 


($2) 


iv / 7? ^ 

7?o KyMDj \yo) 


2/i-3 


1 - 



2^ — 3 


(16) 

where N = {A"^/8TT^)2‘^>^r‘^{y)/{2y - 3). 

Although a full treatment of the semiclassical Ein¬ 
stein equations = —87rG(T^,y + {Tfn,)) lies beyond 
the scope of this work, it is easy to perform a qual¬ 
itative analysis. For small values |^| < 1, the term 
l/(47ra^) f dkk'^lcpkl'^ in Eq. (||) is dominated by horizon- 
size modes, and is small compared to the remaining terms 
in (Too), which are dominated by the larger infrared cut¬ 
off scale rfo {ifo ^ y)- These terms, excluding m^($^), 
are of the form and can be absorbed into the left 

hand side of the (00) Einstein equation leading to: 

3il2~8^G[p™ + im2($2)] (17) 

where G ~ G/(l + 87rG|^|($^)). For ultra-light fields 
1^1 < 0.1 ensures that ($^) grows at a slow rate, satisfying 
constraints on the time variation of G . Consequently, 

8^G($2) at [' (1 + 

\ mpi J \ Tmd J 

(18) 


where ( = fjo/y. (Note that ( = kt,/kmin = exp Hdt, 
where tin is the beginning of inflation, and t* the time 
when a mode entering the horizon today, left the Hub¬ 
ble radius during inflation. In general can be very 
large, log^ 1-) Substituting Hrh/'mpi ~ 10“®, 
Trh ^ lO^^ GeU, Tmd (1 + zmd) x 2.3_x lO-^^ 

1 + Zmd — 23,219 we find that 87rG(4’^) can also 

be very large (since G < G this holds for 87rG($^) as 
well). For large 87rG|^|(4>2) » 1, G ~ l/(87r|^|($2)) and 

Aoff= 8^G(Too) ~ mV2|C| 

= (19) 


We thus find that the energy density of created particles 
behaves like an effective cosmological constant. Further¬ 
more, for ultralight fields m/H ~ 1, can be a sig¬ 
nificant fraction of the total density of the Universe at 
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late times. (Recent observations of high redshift super¬ 
novae suggest Ha > 0, which could be quite large if the 
Universe is flat ll). 

At first glance it may appear strange that the equa¬ 
tion of state of created nonminimal scalars is p ~ —p and 
not the familiar p = p/3 expected for relativistic parti¬ 
cles. This is because the dominant contribution to (Too) 
comes from wavelengths larger than the horizon size, as 
a result (T^jy) — at late times. In a similar 

context it is well known that (Too) for gravitational waves 
(and minimally coupled massless scalars) created during 
inflation is dominated by wavelengths of order the hori¬ 
zon size, because of which the equation of state of relic 
gravitational waves is identical to that of classical matter 
driving the expansion, and in a matter-dominated Uni¬ 
verse, is p = 0 [^. (For wavelengths larger than the 
Hubble radius the distinction between “real” particles 
and vacuum polarization becomes ambiguous.) 

Our analysis indicates that, at a local level, the Uni¬ 
verse may be filled with a quasiclassical, stochastic, 
nearly homogeneous nonminimal scalar field having a 
Gaussian probability distribution with dispersion ($^). 
(This description agrees with the general stochastic ap¬ 
proach to inflationary cosmology 10.) Due to cosmic 
variance, the local value of will generally not equal 
($^), remarkably this does not affect the value of 
since (0) does not depend upon (provided $ is large). 

To summarize, the inflationary production of light non¬ 
minimal particles (^ < 0) leads to an energy density 
which mimics an effective cosmological constant implying 
VLrn < 1 in clustered matter in a critical density Universe. 
Hence, a low density Universe can exist without an “H 
problem” making inflation compatible with observations 
indicating a low value for [0 , 0 . Since nonminimal 
scalars naturally arise in the low energy limit of string 
theory [0, our results may apply to a wider class of 
models than the ones considered here. 

It is a pleasure to thank Yuri Shtanov and Alexei 
Starobinsky for stimulating discussions. 
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